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Abstract 

We measure numerically the distribution of baby universes in the crumpled 
phase of the dynamical triangulation model of 4d quantum gravity. The rele- 
vance of the results to the issue of an exponential bound is discussed. The data 
are consistent with the existence of such a bound. 



One of the more promising approaches to understanding the nature of four- 
dimensional quantum gravity has arisen through models based on summing classes 
of simplicial manifolds - the dynamical triangulation (DT) models |l[ Q . The man- 
ifold is approximated by a set of equilateral simplices whose edge lengths are taken 
to constitute an invariant cutoff. Quantum fluctuations of the geometry are incor- 
porated by constructing a partition function which sums over all possible ways of 
assembling these simplices into a piecewise linear manifold 

Z(k ,k 4 )= e KoN °- K * N *. (1) 

T(S 4 ) 

Here the class of triangulations has been restricted to that of spherical topology. 
The coupling K4 constitutes a bare cosmological constant conjugate to the total 
number of four-simplices (volume) JV4. Similarly, kq plays the role of an inverse 
Newton constant coupled to the total number of zero-simplices (nodes) iVo in the 
triangulation. 

Assuming that we wish to remove the edge length cutoff it is necessary to 
find points in the parameter space of the model at which the mean volume (N4) 
diverges. To see how this may happen consider expanding the grand canonical 
partition function Eq. (1) as a power series in e~ K4 , 

Z(k ,k a ) =J2 n (^ N *) e ~ K4N4 • ( 2 ) 

The coefficients in this expansion are the microcanonical partition functions for the 
system at fixed volume N4. It is these quantities which are estimated in Monte 
Carlo simulations. In two dimensions it is known rigorously that the analogous 
coefficients behave as Q(N2) ~ e K z N2 - that is there is an exponential bound on the 
number of triangulations composed of N2 triangles provided we restrict the global 
topology sufficiently^]. The existence of this bound ensures that the expansion has 
a finite radius of convergence determined by the critical coupling The mean 
volume can then be shown to diverge in power-like fashion as this critical coupling 
is approached. This is the basis for taking the continuum limit. 

In dimensions greater than two the volume dependence of 0(^0,-^4) even 
when restricted to the four sphere is, in principle, unknown. In a previous paper 
we pointed out that the behaviour for small volume is consistent with a super- 
exponential growth Q (N4) ~ g/^iogA^ ||_ This would, at least naively, render a 
continuum limit impossible. Since then two other groups have examined the issue 
on larger lattices and claim strong evidence for an exponential bound @, |||. In light 
of this we have both extended our calculations to larger volumes and in addition 
looked at alternative quantities such as the distribution of baby universes. The 
latter is very sensitive to the volume dependence of Q(kq,N4) and might thus be 
useful in resolving this issue. 

1 In two dimensions this restriction amounts to fixing the genus of the surface. 
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The usual way in which an exponential bound is observed is by looking 
at the volume dependence of the quantity k1(kq,V) = y\ogfl(V, kq) which is a 
numerical estimate for the effective critical cosmological coupling k\ at volume V . 
In Fig. 1 we plot it as a function of the logarithm of the volume to expose any 
logarithmically divergent component to the critical coupling. We show both our 
data (circles) together with the data published in |||, ||. Clearly, these different 
simulations are in agreement within statistical errors. We then fit our data from 
volumes V = 8000 — 128000 using two different functional forms. The straight line 
represents a least square fit to a super-exponential form 

kI(V) = a + l3\ogV , (3) 
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Fit 


a (a') 


0(JT) 


X 2 /d.o.f 


Power 7 = 0.25 


1.242(3) 


-1.23(4) 


2.4/3 


Power 7 = 0.1 


1.389(8) 


-0.68(2) 


0.8/3 


Log 


0.894(8) 


0.025(1) 


3.1/3 



Table 1: The optimal fit parameters as the data for k%(V) is fitted to either super- 
exponential behaviour or a weak power-law convergence. 

while the curve corresponds to a weak power-law convergence to an exponential 
bound of the form 

Kt(V) = a' + ^. (4) 

The fit parameters and quality of the fits are shown in Table 1. Since there 
is so little data we have chosen to do the power fit with two different fixed powers 
7 = 0.25 and 7 = 0.1. Arguments for the former choice are made in Q and 
it corresponds to the curve plotted in Fig. 1. It is clear that both types of fit 
can equally well describe the data. The quality of the fit with 7 = 0.1 appears 
somewhat superior but since the log fit has \ 2 °f order one this should not be taken 
as significant. In || a 7 = 0.25 fit over the same volume range was claimed to 
be substantially better than the logarithm. Our data do not seem to support this 
and we interpret this as simply pointing to the delicacy of deciding between similar 
fits with rather limited data. It is quite possible that many runs at intermediate 
volumes would be useful to resolve this issue. 

Thus while we see that an exponential bound is certainly consistent with the 
existing numerical data at large volumes it is not strongly preferred over the loga- 
rithmic divergence. In light of this we have turned to an analysis of other quantities 
to try to settle the issue. The distribution of baby universes is one such observable 
||. A baby universe is defined as a section of a d-dimensional triangulation con- 
nected to the bulk only through a so-called minimal neck which consists of d + 1 
{d — l)-simplices or faces constituting a boundary of a simplex not already present 
in the triangulation. In four dimensions this is a set of five tetrahedral faces which 
make up the surface of a new simplex and divide the triangulation into two pieces. 
The volume of the baby universe is defined to be the number of simplices in the 
smaller piece. 

The distribution of these baby universes can be computed by considering the 
number of ways a volume V triangulation can be built from a baby of volume B and 
a mother of size V — B by attaching the baby to the mother at some point. This 
gluing operation is effected by identifying one simplex in the baby with another on 
the mother. Thus the distribution takes the form 
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P(B) 



(V - B) Q(kq, V -B) B n{K , B) 

n(K ,v) 



(5) 



Strictly speaking the factors 0,(kq,B) should be replaced with one point functions 
but we shall ignore this unimportant technicality here. The important thing to 
notice is that any exponential factor in £1 cancels out in this formula and P{B) only 
depends on sub-leading corrections - that is it is maximally sensitive to the finite 
volume corrections to coefficients Vl{kq,V). 

With this in mind we have measured the distribution P(B) numerically in the 
crumpled phase of the model when kq = 0. The true volume of the triangulation 
space is most easily estimated here since all triangulations contribute with equal 
weight to the partition sum. Indeed we do not believe it is safe to try to estimate 
the behaviour of Q(kq,V) from simulations at large kq. At such node couplings 
the dominant triangulations correspond to branched polymers whose mean node 
number varies linearly with volume. Such configurations are known to possess an 
exponential bound. The crumpled configurations which predominate at small kq in 
contrast have mean node numbers scaling as some fractional power of the volume. 
At large kq these latter configurations will receive large (as V — > oo) exponential 
suppression relative to the branched polymers from the node term in the action. 

We have simulated the model at four different volumes; 500, 1000, 4000 and 
8000 simplices, using runs of length 10 million sweeps^. We will see that the mea- 
sured distribution falls off exponentially fast which necessitated such high statistics 
runs. This precluded the use of significantly larger lattice volumes in this study. Us- 
ing our previous parameterizations of the finite volume corrections to £l(V) we have 
attempted to fit the data with functional forms corresponding to either logarithmic 
divergence or weak power law convergence 



logP(fl) = a + P((B + 6)\og(B + 5) + (V -B + 5)log(V - B + 5)) , (6) 



The constant 5 is inserted as a phenomenological parameter to reflect sub-leading 
finite size corrections and a and a' reflect an ambiguity in overall normalization. 
In practice we have removed the largest contribution to the latter by dividing the 
measured number of baby universes by the volume V . 

Figure 2 shows the distributions together with a series of curves resulting 
from least-square fits assuming the logarithmic scenario Eq. (6). The fit to the 
largest volume yields a = -2.92(3), = +0.056(1) and 5 = -7(1) with X 2 = 9.6/6 
(per d.o.f.). Fits to the other volumes give consistent results. Notice that we have 
fitted baby universes with size B = 4(n + 1) only (n integer). Baby universes of 
size B = 4 (n — 1) lie on a curve which while yielding consistent fits for (3 is shifted 
by a constant with respect to the first. This effect has been observed before |7j and 

2 One sweep corresponds to V attempted elementary local moves. 



logP(B) 



a 



+ ((P + bf^ + (V-B + 5) 1 - 7 ) . 



(7) 
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Figure 2: log P{B) vs B with a logarithmic fit. 



is presumably the result of finite size effects. We fit only for B > 10 and truncate 
due to poor statistics at B > 50. 

Figure 3 shows the same data now fitted according to the power scenario 
Eq. (7). The best fit in this case yields a' = -0.2(15), /?' = -1.38(5) and 5 = 3(2) 
with x 2 = 6.3/6 assuming 7 = 0.25 as before. At face value then it remains hard 
to differentiate between the two situations. However, notice that the extracted 
value of (3 = 0.056(1) from the log fit is more than twice its estimated value from 
the fits for the effective critical coupling [3 = 0.025(1) (Table 1). In contrast the 
estimate for [)' = —1.38(5) from the power fit is quite close to its value estimated 
earlier (3' = —1.23(4). The relative proximity of the two estimates is particularly 
impressive considering that one is derived from the behaviour of baby universes with 
size less than 8000 simplices while the other is extracted from the critical coupling at 
volumes much greater than 8000. Furthermore, it is clear that the power fit would 
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Figure 3: logP (B) vs B with a power fit assuming 7 = 0.25. 



still hold good if we set 5 = a' = so that such a fit (with a truly minimal number 
of parameters) would do much better than the logarithm. 

Additional information can be obtained by looking at the mean number of 
nodes per unit volume. It is easy to see that this quantity is related to the critical 
coupling through (see eg M) 



Nq\ = 3k% (V,k ) 
V / dK 



(8) 



Thus finite volume corrections to the effective critical coupling result in similar finite 
volume corrections to (Nq/V). In Fig. 4 we show this quantity on a log-log scale 
together with a least-square power fit. While the x 2 °f such a fit is terrible, showing 
that such a simple parameterization is insufficient to describe the data in detail, the 
fit shows that a small power-law correction is again rather well able to account for 
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Figure 4: Mean node number per unit volume with fit assuming a weak power 
convergence. 

the overall structure of the finite volume corrections^]. Notice that any coefficient of 
a logarithmic piece in k% (V, kq) will not contribute since it cannot depend on kq. 

In conclusion, we have presented numerical results which, although not defini- 
tive, are very consistent with the existence of an exponential bound in the dynamical 
triangulation model of 4d quantum gravity. The evidence for this comes both from 
fits to the volume dependence of the critical coupling, an analysis of the baby uni- 
verse distribution in the crumpled phase and the scaling of the mean node number. 
Although individually these quantities are not very conclusive, it is remarkable how 
consistent results are obtained if we assume a weak power convergence. Clearly, it 
is important to strengthen these conclusions both by simulating intermediate lat- 

3 The fit yields 7 = 0.26. 
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tice volumes and perhaps via a high statistics simulation at say volume V = 16000 
directed at probing further into the tail of the baby universe distribution. 

The calculations reported here were supported, in part, by grants NSF PHY- 
9503371, PHY-9200148 and from research funds provided by Syracuse University. 
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